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ABSTRACT 

The photometric evolution of galaxies in a hierarchically clustering universe is inves- 
tigated. The study is based on high resolution numerical simulations which include the 
effects of gas dynamics, shock heating, radiative cooling and a heuristic star formation 
scheme. The outcome of the simulations is convolved with photometric models which 
enables us to predict the appearance of galaxies in the broad band colors U, B, V, R, 
I and K. We demonstrate the effect of the mutual interplay of the hierarchical build- 
up of galaxies, photometric evolution, k-correction, and intervening absorption on the 
appearance of forming disk galaxies at redshift one to three. We also discuss to what 
extend the numerical resolution of current computer simulations is sufficient to make 
quantitative predictions on surface density profiles and color gradients. 
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1. Introduction 

Over the past few years, the advent of a new 
generation of large telescopes and the superb 
imaging quality of the refurbished Hubble-Space- 
Telescope have changed the way we can look at 
the formation and evolution of galaxies. Galax- 
ies are routinely identified at redshift z ~ 1 and 
also objects at redshifts z ~ 3 have been optically 
resolved (Steidel et al. 1996). Theoretical work 
has been guided by the hierarchical clustering hy- 
pothesis. In this scenario, the Cold Dark Matter 
(CDM) model is probably the best known ex- 
ample, structures form as halos of progressively 
larger mass merge, virialize and form new higher 
mass objects. Thus, the galaxy formation pro- 
cess can be modeled from ab-initio conditions 
and model predictions can be rigorously com- 
pared to observational data. Two complemen- 
tary techniques have been shown to be especially 
powerful. 

In so-called semi-analytical studies (Kauffmann, 
White & Guiderdoni 1993; Cole et al. 1994), phe- 
nomenological recipes are used to model the ac- 
cumulation of gas at the center of dark matter 
halos and its transformation into stars. The in- 
volved free parameters are calibrated by match- 
ing present day properties of galaxies, e.g., the 
properties of the Milky Way or the Tully-Fisher 
relation. These techniques are comparably cheap 
as far as computation is concerned and thus al- 
low a large coverage of the parameter space. The 
major strength of the semi-analytical approach 
lies in modeling statistical properties of the whole 
galaxy population as, e.g., the luminosity func- 
tion of galaxies. Its disadvantage is that little 
can be said on details of the galaxy formation 
process or the exact structure and kinematics of 
galaxies. 

Numerical hydrodynamical simulations which 
follow the formation of individual galaxies in de- 
tail fill in that gap. Indeed, advanced numeri- 
cal simulations include all the relevant physical 
processes: gravitation, gas dynamics, adiabatic 



heating and cooling, hydrodynamic shocks, ra- 
diative cooling and star formation (Katz 1992; 
Navarro k White 1993; Steinmetz & Miiller 1994, 
1995; Evrard, Summers & Davis 1994, Tissera, 
Lambas & Abadi, 1997). However, the required 
numerical resolution (~ 1 kpc) together with the 
demanding numerical algorithms and the limited 
computational power of the present generation of 
supercomputers constrain current simulations to 
small box sizes, i.e., each simulation only creates 
very few galaxies. Thus hydrodynamical simula- 
tions can deal with the details of the galaxy for- 
mation process, but can make only limited pre- 
dictions on statistical properties. 

Numerical simulations have also helped iden- 
tifying some of the major problems of the hi- 
erarchical clustering hypothesis: Since cooling 
times scale inversely with density, the dissipa- 
tive collapse of gas should have been more ef- 
ficient at high redshift because the dark matter 
halos present at that time were denser, as was the 
universe as a whole. Thus gas cools and settles 
into rotationally supported disks already when 
the first level of the hierarchy collapses. Conse- 
quently, models predict an overabundance of low 
mass galaxies as well as disk galaxies which for 
a given circular velocity are too massive (White 
& Rees 1978). A related problem is the size 
of present day disk galaxies: Due to the fre- 
quent merging angular momentum is efficiently 
transported from the baryons to the dark-matter 
halo resulting in disk galaxies which are too con- 
centrated compared to present day disk galax- 
ies (Navarro, Frenk & White 1995; Navarro & 
Steinmetz 1997). Feedback processes due to star 
formation, like stellar winds or supernovae, have 
repeatedly been envisioned to solve this over- 
cooling problem (e.g., Dekel & Silk 1986). How- 
ever, a simulation which includes the effect of 
star formation and which is consistent with ob- 
servations has not been performed yet. 

In this paper we merge the approach of numer- 
ical simulations with traditional methods of spec- 
trophotometric synthesis. This hybrid scheme 
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enables us to perform an even more realistic com- 
parison between model predictions and actual 
observations. The structure of the paper is the 
following: In section 2 details of the numerical 
simulations and the star formation scheme are 
provided and the spectrophotometric modeling is 
presented. In section 3 we apply these methods 
to model the appearance of high redshift galaxies 
as they form in CDM like universes. In section 4 
we describe to what extend current simulations 
may even allow us to analyze the detailed struc- 
ture of galaxies at high and low redshift. Section 

5 summarizes the results and draws our conclu- 
sions. 

2. Methods 

2.1. Numerical simulations 

The simulations were performed using Grape- 
SPH (Steinmetz 1996), a hybrid scheme of a 
smoothed particle hydrodynamics (SPH) approach 
with a direct summation N-body routine which 
uses the hardware N-body integrator GRAPE. 
The code is specially designed to follow a three 
component system of gas, stars, and dark matter 
and has high resolution in space and time due to 
the use of individual smoothing lengths and time 
steps. It includes the effects of pressure gradi- 
ents, adiabatic heating and cooling, shocks and 
radiative cooling. 

In the simulations presented here, a star for- 
mation scheme similar to the one in Steinmetz 

6 Midler (1994, 1995, henceforth SMI and SM2; 
see also Katz 1992) is used: Star formation is 
modeled by creating new collisionless "star parti- 
cles" in regions where the gas is locally Jeans un- 
stable and where the cooling timescale is shorter 
than the local dynamical time scale. The orbits 
of these newly formed star particles are subse- 
quently followed in a self-consistent fashion, as- 
suming that they are only affected by gravita- 
tional forces. Young star particles devolve en- 
ergy and metal enriched mass to their surround- 
ing gas, an effect that mimics the energy and 



mass input by supernovae and evolving stars into 
the ISM. The supernova energy is added to the 
thermal energy of the surrounding gas. Input 
in kinetic energy (Navarro & White 1993) is not 
directly implemented but is described by an effi- 
ciency parameter which is of the order of a few 
per cent, an assumption that is justified since we 
preferably concentrate on the formation of lumi- 
nous galaxies like the Milky Way. 

2.2. Initial conditions 

We use idealized initial conditions as pro- 
posed by Katz (1991) and discussed in more 
detail in SM2: We consider an isolated sphere 
on which small scale fluctuations according to a 
CDM power spectrum are superimposed. To in- 
corporate the effects of fluctuations with longer 
wavelengths, the density of the sphere has been 
enhanced and rigid rotation corresponding to a 
spin parameter of A = 0.08 has been added. 
The simulation is identical to that presented in 
SMI and SM2 where a disk galaxy similar to the 
Milky Way forms. These idealized initial condi- 
tions are fairly realistic for the early formation 
epochs (z > 2) when the sphere itself has not yet 
collapsed. The late formation epoch, however, 
is dominated by the collapse of the sphere and 
hence proceeds more smoothly than under more 
realistic conditions. Due to the smooth infall at 
low redshifts the simulation does not suffer from 
the angular momentum problem. Furthermore, 
this model allows for very high numerical reso- 
lution (spatial resolution lkpc, mass resolution 
1O 7 M ) at comparably low computational cost. 

2.3. Simulation outcome 

At any redshift, the simulations provide the 
spatial distribution of star particles. Each star 
particle corresponds to a population of about one 
million stars which have formed in a burst, i.e., 
the whole galaxy is a superposition of several ten 
thousands of mini starbursts. Each of the star 
particles is described by its age (time passed since 
creation) and its metallicity. 



3 



10 



-10 



-30 -20 -10 10 20 30 
X [kpc] 



-40 -20 20 40 
X [kpc] 



Y [kpc] 



[kpc] 



Fig. 1. — Upper left: distribution of star particles 
at redshift z = 3.5. The sidelength of the plots 
corresponds to 100 pixels of the HST/WFPC2. 
Upper right: face-on view of the galaxy at red- 
shift z = 1. Lower left: edge-on view of the 
galaxy at redshift z = 1. Lower right: age of star 
particles at z = 1 as a function of their distance 
to the galaxy center. 

The actual distribution of star particles is 
shown in figure |]. Details on the formation his- 
tory of the model galaxy of its kinematical prop- 
erties as well as of its chemical evolution can be 
found in SMI and SM2. Here we concentrate on 
some features important for the discussion later 
in this paper. At redshift 3.5, two protogalac- 
tic clumps can be seen which are going to merge 
within the following ~ 10 8 years. At lower red- 
shifts, the stars originating in these two clumps 
can be found in the bulge of the galaxy. At z = 1, 
the model galaxy already resembles a Milky Way 
type galaxy with a bulge-to-disk ratio of 1:3 in 
mass (SM2). In the lower right panel of figure 
[l], we plot the age of the star particles as a func- 
tion of the distance to the galaxy center. The 
central 1-2 kpc are dominated by old stars with 



an age between 3 and 4 Gyr. The majority of 
these star particles has a metallicity of above so- 
lar (SMI). There is little offset between the age 
of the oldest stars (3.9 Gyr) and the age of the 
universe at that redshift (4.6 Gyr). This feature 
is common to any hierarchically clustering uni- 
verse in which the first stars form in small clumps 
which collapse at high redshifts. The actual value 
here is only an upper limit since (i) the limited 
numerical resolution cannot follow star forma- 
tion in clumps with less then ~ 10 6 M Q and (ii) 
the amount of substructure is underrepresented 
due to the idealized initial conditions. At larger 
radii (r > 2 kpc), younger stars with an age of 
about 1 Gyr and below are dominating which re- 
flects ongoing star formation in the gaseous disk. 
Due to the vacuum boundary conditions, the evo- 
lution of the galaxy is fairly quiescent at red- 
shifts below one. Gas in the disk is successively 
transformed into stars. Due to the still insuffi- 
cient numerical resolution, however, encounters 
between disk stars and dark matter particles act 
as a kinematic heating term. While at z = 1 
the velocity dispersion of disk stars is still small 
(« 10 km/sec), the stellar disk is successively 
heated up between z = 1 and z = up to ve- 
locity dispersions near 100 km/sec at z = 0. We 
therefore focus our analysis mainly on redshift 
one, when the amount of artificial heating is still 
negligible. 

2.4. Spectral synthesis 

As mentioned above, each star particle, which 
has a typical mass of a few million solar masses, 
corresponds to a population of several million 
stars which have been formed at the same time 
tf from a gas cloud with a metallicity of Z. A 
star particle can thus be considered as a pop- 
ulation of stars which have formed in a burst- 
like manner. The luminosity evolution of such a 
burst is then followed by an evolutionary spec- 
tral synthesis model (Fritze-von Alvensleben &; 
Gerhard 1995; Fritze-von Alvensleben & Burk- 
ert 1995; Einsel et al. 1995): The evolution of 
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this instantaneously formed population of stars 
is followed in the Hertzsprung-Russell diagram 
(HRD) using results of stellar evolution theory. 
Sets of stellar tracks (Maeder 1990) for five differ- 
ent metallicities (Z = 1(T 4 , 10~ 3 , 4x 10~ 4 , 1(T 2 , 
4 x 10 -2 ) and for 27 stellar masses in the range of 
0.15 M© to 60 M© have been compiled by Einsel 
et al. (1995). A Miller-Scalo (1979) initial mass 
function (IMF) is assumed with upper and lower 
mass limits of 60 M© and 0.15 M©, respectively. 
The spectrum of an individual star burst is syn- 
thesized from Kurucz's (1992) library of stellar 
spectra, luminosities in various filter bands are 
obtained by convolution with filter response func- 
tions (see also Leitherer et al. 1996, Fritze-von 
Alvensleben & Kurth 1996). 

The discrete set of stellar masses for which 
stellar evolution tracks are available together 
with the burst-like star formation history can 
cause some problems for the late time evolution 
of the population. In the original implementa- 
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tion, all stars within a mass interval 
have been assigned to the temperature and lu- 
minosity of the stellar track of a star of mass 
mj. For masses below 2M©, the lifetime on the 
red giant branch (a few 10 8 yr) is, however, short 
compared to the difference in the lifetime of stars 
on neighboring evolutionary tracks (10 9 yr). This 
causes all stars within a given mass interval to 
synchronously enter the red giant branch. Since 
the luminosity increases hundredfold, and since 
every mass interval contains a substantial frac- 
tion of the total mass of the population (a few per 
cent), pronounced spikes in the luminosity can be 
seen (figure ^, dotted line). Such a discreteness 
dominated behavior is unknown for extended star 
formation histories (as e.g. in models of spiral 
galaxies), since these spikes are convolved (and 
thus smoothed) with the star formation history. 

These discreteness effects in the luminosity 
evolution of burst-like formed populations are 
circumvented by methods which are based on 
isochrones (see, e.g., Bruzual & Chariot 1993). 
We propose an alternative method that avoids 
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Fig. 2. — Luminosity evolution of a stellar popu- 
lation formed within 2.5 x 10 7 years with metallic- 
ity Z = 1 X 10 -4 . The dotted line shows the evo- 
lution as computed with the original code. The 
luminosity indicated by the dashed line is calcu- 
lated with the averaged stellar evolution tracks 
and smoothed subsequently (solid line). The 
numbers right of the peaks assign the mass rrii in 
M©of the star that reaches the red giant branch 
at that time. 

the construction of isochrones. The luminosity 
evolution of an individual star is split into two 
parts, main sequence and post main sequence. 
As can be seen in figure ||, the post main sequence 
evolution between neighboring masses is almost 
identical, concerning the time dependence as well 
as the luminosity evolution. The main sequence 
evolution Zms (t, m ) has also a very similar shape, 
but the evolution is faster for the higher masses 
and the star is also brighter. These two effects 
can, however, be compensated by the following 
ansatz for the luminosity evolution of a star with 
mass m + Am: 

m + Am 



Ims (t, m + Am) = l us (t', m) 

\ m 

(1) 

In this ansatz, a is the exponent of the mass- 
luminosity relation, typical values being a » 3 
for masses less than a few M© . The rescaled time 
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t' can be inferred by rescaling t with the ratio of 
the life times on the main sequence tms> i-e., 
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with j3 ~ — 2; stars with a higher mass reach a 
given evolutionary stage faster. For each mass to 



within the mass interval 
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minosity evolution can thus be calculated. The 
local value of the slope of the mass-luminosity 
and mass-lifetime relation, a and /3, can be 
drawn from neighboring masses mj. In figure || 
we demonstrate the success of our interpolation 
scheme by comparing the evolution of a star of 
mass M = 1.6 Mq with the evolution as inferred 
from a star with M = 1.4 M and M = 1.8 M 
using the interpolation scheme above. Differ- 
ences in the main sequence evolution are almost 
undetectable and also the post main sequence 
evolution differs only by a very small time shift 
(less than 10 per cent of the lifetime on the post- 
main-sequence) . 
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Fig. 3. — Comparison of the luminosity evolu- 
tion of three stars of mass M = 1.4 Mq, 1.6 M© 
and 1.8 Mq, respectively. The dashed line has 
been obtained by interpolating between the M = 
1.4 M & and 1.8 Mq track as described in the text. 

By averaging these approximated stellar evo- 
lution tracks in each mass interval and by adding 



up the luminosities according to the initial mass 
function (IMF) the total luminosity of the stel- 
lar population is obtained as a function of time. 
The luminosity of populations which formed in a 
single star burst now evolves much smoother (see 
figure [2|, dashed line) . The remaining small wig- 
gles result from discontinuities in the evolution 
near the interface of two adjacent mass intervals. 
We remove these wiggles by averaging over time 
intervals of 1.5 x 10 8 years, as indicated by the 
solid line in figure |2[ 

The luminosity evolution in the U and K band 
of a single star burst is shown for 5 different 
metallicities in figure |I[ Data for arbitrary metal- 
licities Z can then be obtained by interpolation 
in log Z. 

2.5. k correction 

For apparent magnitudes and colors of galax- 
ies at high redshifts the modifications due to the 
redshift of the spectrum have to be taken into 
account. 

The cosmological redshift causes the spectral 
energy distribution Fq (A) of a galaxy in its rest- 
frame at redshift z r to be detected in the ob- 
servers restframe as F (A) with: 



F (A, Zr) 



F (A/(l + Zr)) 
(1+Zr) 



(3) 



The apparent magnitude in some a X" filter 
band is thus 



m x = B DM (z r )-2. 5 log 



' VO VM X 0+k(z r ) 



L 



x©/ 



(4) 

with the luminosity Lxo in the galaxy's rest 
frame, the bolometric distance modulus 

BDM (z r ) = 5 log {D L («,.)/Mpc) + 25, 

and the luminosity distance (Mattig 1958) 

D L (Zr) 



"77 — 2 { Z rV0 + (<?0 - 1) 
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1 + J{2q z r + 1) 



(5) 
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Fig. 4. — Luminosity evolution in the U and K 
band of a single age stellar population for the 
metallicities Z = 1 x 10 -4 , 1 x 10 _3 ,4 x 1(T 3 , 1 x 
1CT 2 , and 4 x 1CT 2 . The total mass of the pop- 
ulation is 1 x 1O 7 M . 

The k correction is calculated from the spectral 
energy distribution F (A) and the filter response 
function Sx (A) (Lamia 1982) according to the 
following equation (Oke & Sandage 1968) 



k 



2.51og(Lxobs./£xo) (6) 
l + z r )f™F (X)S x (X) dX 



2.5 log 



(7) 



f™F (X/(l + z r ))S x (X) dX 

2.6. Creating artificial images 

From the luminosity of the set of star particles, 
artificial images can be easily generated. The lu- 



minosity in a given pass band filter can be con- 
verted into number of photons emitted per unit 
time per solid angle. This number is corrected by 
intervening absorption due to neutral hydrogen 
using the models of Madau (1995). By specify- 
ing exposure time and pixel size, the number of 
photons received by each pixel is calculated. The 
pixel map then is convolved with the numerical 
resolution of the simulation and the sky back- 
ground is added. Finally it is convolved with the 
point-spread function of the instrument. The re- 
sulting array contains the expectation value for 
each pixel to receive a photon. This number is 
then random-sampled. Together with the detec- 
tor characteristics (quantum efficiency, readout 
noise, dark current, gain) an artificial image can 
be created which subsequently can be analyzed 
in the same way as raw observational data. 





Fig. 5. — Artificial HST images of a forming 
galaxy at z = 3.5. Images are 50 pixels across, 
"exposure time" is 30000 sec. Top row: [/-band; 
bottom row: /-band. Left column: rest frame; 
middle column: including k-correction; right col- 
umn: including k-correction and absorption due 
to intervening hydrogen. 

3. Applications 

3.1. Appearance of high redshift galaxies 

In this paragraph we discuss the appearance 
of the simulated galaxies in an exposure as it 
is typical for HST/WFPC2 observations. De- 
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Fig. 6. — Artificial HST images of a forming 
galaxy at z = 1. Images are 50 pixels across, 
"exposure time" is 30000 sec. Top row: U -band; 
bottom row: /-band. Left column: face on, 
rest frame; middle column: face on, including 
k-correction; right column: edge-on, including k- 
correction 

tector characteristics and sky background have 
been taken from the WFPC2 handbook and a 
generous exposure time of 30000 s has been speci- 
fied. For the adopted cosmological model (0 = 1, 
Hq = 50 kms -1 Mpc" 1 ) and redshift range (z = 
1 — 4), the formal numerical resolution of the 
simulation surpasses the pixel resolution of the 
WFPC2 camera by about 50%. 

In figure [| and ^ we show a U and /-band 
image of the simulated galaxy at redshifts z = 
1 and z = 3.5. Each image has a side length 
of 100 pixels, corresponding to 70kpc (z = 3.5) 
and 82kpc (z = 1), respectively. Grayscales have 
been chosen logarithmically to cover the interval 
between a signal-to- noise of 1 to 30 (U band) and 
100 (J band). 

In the restframe [/-band image of the z=3.5 
galaxy, only the most massive clump is visible. 
It has an apparent magnitude of mjj = 24.8. Ap- 
plying the k-correction makes the clump barely 
visible {mjj = 25.7). At these redshifts, the 
[/-band is centered at wavelengths shortward 
of 912 A in the galaxy's restframe where stel- 
lar emission is depleted due to absorption in 
the stellar atmosphere. Absorption due to inter- 
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Fig. 7. — Evolution of the model galaxy in the 
U—B vs. V — R plane. Solid line: model galaxy, 
including k-correction and intervening absorp- 
tion; dashed line: model galaxy, k-corrected, 
no absorption; dotted line: E galaxy; dashed- 
dotted: Sd galaxy (from Moller et al. 1997). 

vening neutral hydrogen suppresses the [/-band 
emission of the proto galaxy by further 3.5 mag 
(mjj = 29.2) and makes it invisible even in long 
exposures. 

In the /-band, the k-correction has the oppo- 
site effect: in its restframe, the apparent magni- 
tude is mi = 24.6 and also the less luminous pro- 
genitor is barely visible. The luminosity of the 
galaxy is dominated by stars which have formed 
within the last few 10 s yr and also massive stars 
M > 8M are abundant. Thus, these lumps 
strongly emit in the near UV (~ 4000 A). The 
cosmic redshift translates these frequencies to the 
/-band resulting in a negative k-correction of 1.7 
magnitudes (m/ = 22.9). The galaxy thus has 
a brightness compatible to the brightest galax- 
ies in the Steidel et al. (1996) sample. Since at 
z = 3.5 the /-band corresponds to wavelengths 
longward of 1217 A, the luminosity is not affected 
by intervening absorption. 

The restframe [/-band image of the z = 1 
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galaxy consists of a luminous bulge of stars with 
a typical age of 3 — 4 Gyrs (see figure [j], lower 
right panel) and some light from star forming re- 
gions in the disk. At redshift 1, however, the 
[/-band corresponds to restframe wavelengths 
near 1800 A where the emission of old stars is 
small. Thus the bulge component virtually dis- 
appears after applying the k-correction. Simi- 
lar to the z = 3.5 /-band, k-correction enhances 
the star forming regions in the disk. Young 
massive stars emit strongly at restframe wave- 
lengths near 2000 A which are redshifted into the 
[/-band. Similar effects can be seen in the I- 
band image: Since the /-band corresponds to a 
restframe wavelength of 4000 A, also the light of 
the old bulge component is clearly visible, al- 
though slightly dimmed. For the disk popula- 
tion the k-correction favors shorter wavelengths 
at which the contribution of newly formed stars 
is larger. Consequently, the k-corrected /-band 
image looks more structured and the regions of 
active star formation are more pronounced. 

In figure the evolution of the model galaxy in 
the U — B vs. B — R plane is shown and, for com- 
parison, the evolution of a Sa and Sd galaxy us- 
ing traditional spectrophotometric models (see, 
e.g., Bruzual & Chariot, 1993, Fritze-von Al- 
vensleben, 1993) is plotted. 

3.2. Surface brightness, color gradients 

In figure || we show surface brightness profiles 
and color gradients. The radial profile clearly 
shows an exponential behavior for r > 5 kpc and 
a steepening (bulge) at smaller radii. The non- 
exponential bulge component is readily visible in 
the z = plot since the surface brightness of 
the disk has dimmed by about one magnitude. 
A comparison of the k-corrected (figure 
galaxy at z = 1 with its restframe (figure | 
die) exhibits a substantial (slight) enhancement 
of the K (I) luminosity (i.e., the k-correction 
is negative), while the V and U luminosity is 
dimmed. 

At z = 1 the galaxy exhibits very strong red- 
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Fig. 8. — Left column: Surface density versus 
radius. The solid line corresponds to the bolo- 
metric luminosity, the dashed, dotted, dashed- 
dotted and dashed-triple-dotted line correspond 
to the broadband colors U, V, I and K, respec- 
tively. Right column: color indices U— V (solid), 
B-V(dotted), V- /(dashed) and V - K (dashed- 
dotted) versus radius. From top to bottom: 
z = 1 including k-correction, z = 1 restframe 
and z = 0. 

dening within the inner 5 kpc, about 2 magni- 
tudes inV — K and 1 magnitude in the other col- 
ors. Except U — V, the drop is much weaker in 
the restframe. This reddening is consistent with 
the dimming of the bulge component of the z = 1 
galaxy due to k-correction as described in section 
3.1. The bulge is dominated by old stars (see 
figure H|, lower right panel) emitting only little 
near 1800 A which corresponds to the U band at 
z ~ 1. The k-correction enhances the K and / 
luminosity while decreasing the U and V lumi- 
nosities and thus results in the strong reddening 
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as can be seen in figure || Only a slight color 
gradient is visible at radii r > 10 kpc in the rest- 
frame as well as in the k-corrected profile. Also 
the effect of the k-correction is much more mod- 
erate, between mag in B — V, 0.6 mag in U—V 
and 1.5 mag in V—K. At these radii, stars are 
still being formed resulting in a much flatter spec- 
tral energy distribution and thus a much weaker 
k-correction. 

At z = 0, only very weak color gradients are 
visible. The bulge differs from the disk by about 
0.1 mag, only in U — V a slightly stronger color 
gradient is visible. The change of color is cer- 
tainly much weaker than the observed one of 
about 0.2 mag in B — V and 0.5 mag in U—V 
as in galaxies like, e.g., Andromeda (Segalovitz 
1975, Josey & Arimoto 1982). One reason for this 
discrepancy is the neglect of dust absorption. As 
we will show in the next paragraph, however, the 
lack of significant color gradients may also be an 
artifact of our star formation scheme and of the 
vacuum boundary conditions. 

In figure |9] we investigate to what extend the 
color gradients at z = 1 and z = reflect age 
and metallicity gradients of the stellar disk. First 
we estimate the statistical error of the color in- 
dices. For that purpose we perform a bootstrap 
analysis, i.e., we draw 100 random samples from 
the distribution of stars and determine color as a 
function of radius. We use the [10; 90] percentile 
values as the lower and upper error limits, i.e., 
only 10% of our bootstrapped samples have a 
color bluer than our lower error limit and only 
10% of our bootstrapped samples have a color 
redwards of our upper limit. As can be seen, col- 
ors can be measured with a statistical error of 
about ±0.02 to 0.05 mag at z = 0. At z = 1, 
errors are at the level of about ±0.1 to ±0.2. 
These errors strongly depend on whether or not 
stars are actively formed. While colors are fairly 
accurate for passively evolving populations, er- 
rors are much larger for regions including young 
stars. This can be easily seen in the z = U — V 
color at a radius of 12 kpc. This radial bin con- 
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Fig. 9.— k-corrected V — K (left) and U-V 
(right) as a function of radius for the model 
galaxy at redshift z = 1 (top) and z = (bot- 
tom). Stars correspond to the model galaxy, er- 
ror bars are [10; 90] percentile values based on a 
bootstrap analysis (for details, see text). The 
dashed lines correspond to a model with con- 
stant metallicity while stellar ages are taken from 
the simulation. The dotted lines correspond to a 
model with constant stellar ages while metallici- 
ties are taken from the simulation. A comparison 
of the actual values with the dashed and dotted 
lines thus demonstrates the effect of metallicity 
and age on the color gradient. 

tains very few very young stars. The color thus 
strongly depends on how many of these young 
star particles happen to be included in the boot- 
strapped sample. 

We also compare the radial color profiles with 
those obtained assuming a constant metallicity 
for the star particles. For the z = 1 model the 
color near the center is consistent with a metallic- 
ity of Z ~ Zq, while the color at 15 kpc requires 
a metallicity much below 0.5 Zq. No constant 
metallicity model fits the entire color gradient. 
The metallicity gradient inferred from the color 
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is consistent with the directly measured metallic- 
ity gradient of 0.05 kpc -1 . The color gradient is 
not consistent with the assumption of a constant 
age. U — V near the center requires a stellar age 
of 3Gyr, while at radii of 15 kpc an average age 
of a few 10 8 years seems to provide the best fit, 
being consistent with the age distribution shown 
in figure |]. A comparison of U—V and V — K 
also demonstrates that the age distribution at a 
particular radius is inconsistent with the assump- 
tion of a constant age: For example, at radii be- 
low r = 5 kpc the age which fits V — K is more 
than 2 times larger than the best fitting age for 
U-V. 

At z = 0, this picture changes. The color 
gradient is consistent with a constant metallicity 
of about 0.5 Z & and only a very slight age gra- 
dient is visible. Z = 0.5 — 1.0 Z & is the average 
metallicity of Milky Way bulge stars (Mc William 
& Rich 199), Z = 0.5 Z & compares well with 
the stellar metallicities in Sb models of Moller et 
al. (1997). However, the actual metallicity gradi- 
ent in this model is much weaker (< 0.01 kpc -1 ). 
A comparison of the z = 1 and z = model 
thus provides a further argument for the lack of 
a significant color gradient at z = and the too 
small metallicity gradient. As shown in SMI, the 
model galaxy also exhibits a too large fraction of 
low metallicity stars (the so-called G-dwarf prob- 
lem). The reason for the failure to reproduce a 
significant metallicity gradient can be found ei- 
ther in the chosen star formation law or in the 
vacuum boundary conditions or both. Between 
z = 1 and z = the disk is dynamically inac- 
tive. Only gas is continuously transformed into 
stars. Thus each fluid element can be consid- 
ered as a closed box model. The metallicity of 
a closed box model increases monotonically until 
the yield is reached. The yield is, however, iden- 
tical for each fluid element since the same IMF 
has been used. The almost constant metallicity 
at z = thus reflects the fact that the metallic- 
ity of each fluid element has already reached its 
asymptotic value. This artifact can be avoided 



by: (i) decreasing the star formation efficiency; 
in that case a metallicity gradient forms since the 
metallicity of low density gas at larger radii does 
not yet reach its asymptotic value; (ii) inflow of 
gas at redshift z < 1, i.e., non- vacuum boundary 
conditions. We argue that (ii) will be a neces- 
sary entity, not only since such a late infall is a 
generic prediction of any hierarchically clustering 
scenario, but also because it may help to cure the 
G-dwarf problem. 

4. Discussion and conclusions 

We have simulated the formation of a Milky- 
Way type galaxy in a hierarchically clustering 
scenario using idealized initial and boundary con- 
ditions. The simulation outcome has been an- 
alyzed with spectrophotometric methods. The 
most important results are: 

1. A hybrid scheme of high resolution numer- 
ical simulations of galaxy formation and 
spectrophotometric modeling is presented. 
We demonstrated the capabilities of such 
an approach. In particular, we were able 
to model galaxies with an appearance qual- 
itatively similar to observed galaxies at red- 
shifts between z = and z = 4. Using the 
appearance of z = 1 and z = 3.5 galaxies 
we also could nicely demonstrate the mu- 
tual interplay between the galaxy forma- 
tion process, k-correction and intervening 
absorption. 

2. The numerical resolution of present day 
computer simulations enables us to study 
the detailed structure of galaxies func- 
tion of redshift, such as surface brightness 
profiles and color gradients. Numerical res- 
olution is high enough and statistical errors 
are sufficienly small to infer whether a sig- 
nificant metallicity or age gradient of the 
stellar population are present. 

3. The properties of a model galaxy are among 
others dependent on the particular star 
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formation scheme, IMF, supernovae feed- 
back. We are, however, optimistic, that a 
rigorous comparison of high-resolution nu- 
merically simulated galaxies with the vast 
amount of detailed observations at differ- 
ent redshifts may allow us to constrain 
and tighten the parameter space. Such 
an ansatz is very similar to the calibration 
step performed in phenomenological mod- 
els. The work presented in this paper is 
only a very first, but necessary step into 
that direction. 

4. Although this work is based on a single 
numerical simulation and a particular star 
formation model, it already provides some 
interesting physical results. The models 
fail to produce a significant color gradient 
in spiral galaxies at low redshifts and are 
subject to a G-dwarf problem. Both prob- 
lems arise from the idealized boundary con- 
ditions, in particular the lack of any late in- 
fall of gas. Other explanations include, but 
are not limited to, the neglect of reddening 
due to dust. 
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